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The top quark cross section close to threshold in e'^ e~ annihilation is computed including the 
summation of logarithms of the velocity at next-to-next-to- leading- logarithmic order in QCD. The 
remaining theoretical uncertainty in the normalization of the total cross section is at the few percent 
level, an order of magnitude smaller than in previous next-to-next-to-leading order calculations. This 
uncertainty is smaller than the effects of a light standard model Higgs boson. 



Top quark pair production close to threshold is a ma- 
jor part of the top quark physics program at a future 
lepton collider. Near threshold the top quark velocity v 
is small and the presence of Coulomb singularities make 
the summation of terms proportional to ag/v mandatory. 
The top quark width, Tf ~ 1.5 GeV 3> Aqcd, serves as 
an infrared cutoff, allowing for the use of perturbative 
methods to calculate the non-relativistic top-antitop dy- 
namics to a high degree of precision. 

Recent NNLO QCD calculations of the total cross sec- 
tion Gti showed that the top quark mass can be deter- 
mined with uncertainties below 200 MeV, but indicated 
that the strong coupling, the top Yukawa coupling, and 
Ft cannot be measured with good precision due to large 
theoretical normalization uncertainties of about 20% [Q . 
A common feature of all recent NNLO QCD calculations 
is that they are fixed order calculations, i.e. the running 
from the hard scale down to the non-relativistic scales 
which govern the dynamics of the top-antitop system was 
not taken into account. For tt, v ^ 0.15 so that loga- 
rithms of ratios of mt ^ 175 GeV, rritv ~ 25 GeV and 
rritV^ ^ 4 GeV are not small, and the renormalization 
group evolution is significant. 

In this letter we calculate the photon induced tt pro- 
duction cross-section in the framework of vNRQCD using 
the velocity renormalization group (VRG) Q| . The calcu- 
lation includes a summation of logarithms of the ratios of 
the scales rrit, mtv and rritv'^ at next-to-next-to-leading 
order. The VRG-improved computation has a scale un- 
certainty of 2-3% at the peak of the cross section (with- 
out initial state radiation and beam smearing effects). 
An improvement in the convergence of the expansion is 
also found. Measurements of the strong coupling, the 
top Yukawa coupling, and the top quark width appear 
feasible with small theoretical uncertainties. 



The expansion for the normalized cross section R 



takes the form 



R = V (") ^'^^ Inf)* 




(LL) 
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The free quark cross-section is of order v. The Coulomb 
summation of powers of a^/v, and the VRG summation 
of powers of aslnv, are the sums over k and i, respec- 
tively. Terms in the cross section at leading-logarithmic 
(LL), next-to- leading- logarithmic (NLL), and next-to- 
next-to-leading-logarithmic (NNLL) order are indicated 
in Eq. (0). 

vNRQCD D is an effective field theory which describes 
non-relativistic heavy quarks with mass m interacting 
with soft gluons with four-momenta k^^ ~ mv and ul- 
trasoft gluons with k^^ ^ mv^ , where mv^ is larger than 
Aqcd- For loop integrations over soft energies and mo- 
menta the MS subtraction scale is fis — mv, whereas for 
loop integrations with ultrasoft energies and momenta it 
is /i(7 = niiP . The subtraction velocity v correlates /is 
and /iy and is used instead of the MS scale parameter /z. 
The correlation is mandatory since energy and momen- 
tum are related through the quark equations of motion. 
In QED this has been shown to be necessary to repro- 
duce through running the (Ina)*^, fc > 2 contributions in 
Lamb shifts, hyperfine splittings and corrections to the 
ortho and para positronium decay widths Q . 

Lowering v from v — 1 to the quark velocity v, sums 
all large logarithms involving the soft and ultrasoft scales 
into the Wilson coefficients of the vNRQCD operators. 
At NNLL this includes logarithms originating from ra- 
diation effects. Once v is lowered to v, power counting 
shows that matrix elements with ultrasoft gluons do not 
have to taken into account at NNLL for the description 
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of a heavy quark pair in a color singlet state. Thus, af- 
ter lowering v to order ag (since ?; ~ in a Coulombic 
system) the NNLL equation of motion of a color sin- 
glet quark-antiquark system is a conventional two-body 
Schrodinger equation. 

In momentum space the NNLL Schrodinger equation 
reads 



,2 p4 



E 



= (27r)"<5(")(p-p'), 



(2) 



where m is the heavy quark pole mass, E = (y^ — 2m), 



3 - 2e, and 



-(47r)-^d"q/(27r)". The po- 



tential for the quark-antiquark pair in a Si state (rele- 
vant for top production through a virtual photon) is 



where (og - 



Vcip, q) + 14 (P, q) + Vs{p, q) -f V;(p, q) , (3) 
= as(Ms), L = ln(fcVA*|), k = p q) 
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The potentials arise from the four-quark matrix elements 
of potential-type operators and from time-ordered prod- 
ucts of operators describing interactions with soft glu- 
ons. The coefficients of the potentials at the hard scale, 
Vi{iy = 1), are obtained with on-shell matching. For 
v V the coefficients contain the summation of all Inw 
terms. The velocity counting of each potential in Eq. ^ 
is equivalent to the counting in the Schrodinger equation 
in previous NNLO calculations (see e.g. Ref. jj]). The 
coefficient Vc{i^) in the Coulomb potential, Vc, was de- 
termined in Ref. |^ at NNLL order. The second term in 
Eq. (^ contains the one and two-loop corrections to the 
Coulomb potential In vNRQCD it arises from the 
time-ordered product of the lowest order operators de- 
scribing the interaction of quarks with soft gluons. The 
couplings in this time ordered product are as (/is) and 
evolve with the QCD /3-function. The potential Vk leads 
to terms in the cross section that are w^-suppressed. The 
coefficient has a matching value of order , so its 
NLL evolution from Ref. |^ is needed. The potentials 
and Vr also lead to terms in the cross section that are 
f ^-suppressed and have order as coefficients V2.s,r gener- 
ated at tree level. Their evolution is only needed at LL 
order and can be found in Ref. M . 



To describe vector current induced quark-antiquark 
production close to threshold at NNLL order we also need 
the Wilson coefficients of the dimension-three ^Si cur- 
rents at NNLL order and the corresponding dimension- 
five currents at LL order. The non-relativistic current for 
production is Jp = ciOp^i + C20p.2, where 



Op i = ijjp'' cr {{(72) X-p , 

Op 2 = ^ ipp^ {p^ t) {ia2) X- 



(8) 
(9) 



Spin and color indices are suppressed. There is an- 
other dimension-five current describing D-wave produc- 
tion which does not contribute at this order. The cor- 
responding annihilation currents Op, 1,2^ are obtained by 
complex conjugation. The matching condition for the 
Wilson coefficient ci at the hard scale needs to be known 
at order a^, and the Wilson coefficient C2 needs to be 
known at the Born level. The value of ci{v = 1) can 
be determined from matching the two-loop result for the 
quark-antiquark production amplitude close to threshold 
in full QCD [|l^ to the corresponding amplitude in vN- 
RQCD. We find 



ci(l) = l-2Cf^^ + a2(TO) 



In 2 25 2 
T2" ^ 24 ^ ^ 



^CaCf[ ln2- 1 



(10) 
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The two- loop result for ci(l) is scheme-dependent and 
our result differs from the hard contribution obtained 
from the threshold expansion [p^p|. 

The LL anomalous dimension for ci is zero. The evo- 
lution of ci for 1/ < 1 at NLL order has been determined 
analytically in Ref. |H| by solving H 
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The resummed logarithms in ci{i' = Us) at NLL order 
include the sizeable negative normalization corrections 
cx a^ln^Qfs found in Ref. At NNLL order we find 

that there are no additional contributions to the anoma- 
lous dimension for ci, so Eq. ( p^ remains valid. However, 
solving for the full NNLL ci{v) requires the NLL values 
of Vc, V2, Vr and Vs, and the NNLL value of Vfe. Besides 
Vc, the coefficients at this order are not yet known. For 
our NNLL solution we will use ci(l) at 0{a1) with NLL 
evolution in v. From previous experience in weak decays, 
it is not expected that the NLL (NNLL) results for V2, Vr 
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and Vs (Vfe) will deviate significantly from the LL (NLL) 
ones. Thus, our result for the top-antitop cross section 
should yield a realistic estimate of the theoretical uncer- 
tainties at NNLL order. The LL result for C2 is generated 
by ultrasoft gluons renormalizing Op i and reads 



32, /as(mi/2) 
as{m) 



(13) 



The first logarithm in this series agrees with the loga- 
rithm in the matching calculation in Ref. [ p^ . 

In full QCD the expression for the normalized photon- 
induced cross section for heavy top-antitop production 
R = Oti/'y^+fi- at the cm. energy ^/s is 



167r 
R = Im 

9s 



-I d^xe''>-^{0\Tr{x)jMG) 



, (14) 



where q = {^/s, 0) and is the vector current that 
produces a top-antitop pair. In vNRQCD at NNLL or- 
der the vector current correlator is replaced by the cor- 
relators of the non-relativistic currents Op^i and Op, 2 
evaluated for ly ^ a^. The correlator of two Op i cur- 
rents is proportional to the coordinate space Green func- 
tion G(0, 0) obtained from the NNLL Schrodinger equa- 
tion (Q) and the correlator of Op i and Op 2 is propor- 
tional to {E/m)G{0,0). To determine G(0,'o) at NNLL 
order we use a combination of numerical and analytic cal- 
culations [|^. The ultraviolet-finite contributions from 
the Coulomb potential (^) (denoted by in Eq. ( [l5| ) be- 
low) are determined exactly using numerical techniques 
developed in Ref. [|5|. The ultraviolet-divergent contri- 
butions from V5, Vr, Vk and the kinetic energy correction 
(denoted by SG^ , SG^, 5G^ and JG"^'", respectively) are 
computed in perturbation theory with dimensional regu- 
larization in the MS scheme. This is required for consis- 
tency with the scheme used to compute the matching and 
running of the Wilson coefficients. By power counting, 
only a single insertion of these potentials is included. In 
deriving these results we have included the counter terms 
generated by renormalizing the Op i current These 
counter term graphs are sufficient to cancel all subdi- 
vergences. The remaining overall divergences are of the 
form 1/e and u^/e and are cancelled by renormalizing the 
time ordered product of currents. The final result for the 
NNLL cross section in vNRQCD is 



R 



NNLL 



^ Im| cl\{\- w^) G^ + (V2 + 2Vs) SG^ 
L L 



+Vr 5G' + Vk SG'' + SG 



kin 



+ 2ciC2V^G^^, (15) 



where v is the top quark velocity and the dependence 
of the Wilson coefficients on 1/ is suppressed. The top 
quark width is implemented by shifting the energy into 
the positive complex energy plane by iVt JTot , which leads 
to the following expression for the top quark velocity 
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FIG. 1. The upper panel shows the normalized pho- 
ton-induced cross section for ti at LL order (dotted), NLL 
order (dashed), and NNLL order (solid) in the IS mass 
scheme with v = 0.15, 0.2, and 0.4. The lower panel 
shows the v dependence of the contributions in Eq. (|l^) to 
^NNLL^yj = 350 GeV). The contribution is shown by 
the dotted line, the sum of SG^-''-^'" terms by the dot-dashed 
line, the terms involving SC^ by the dashed line, and the sum 
by the solid line. 



G — 2mt + iVt 
nit 



(16) 



We emphasize that Eq. ([ij) does not provide a consis- 
tent treatment of the top quark width beyond next-to- 
leading order. This can be seen from the presence of in- 
dependent terms proportional to asFt/TOj ln(— w/j/) in 
Eq. ( p^ ) which are parametrically of NNLL order [Q. 
Conceptually, they indicate that further renormalization 
procedures are required in a consistent treatment of elec- 
troweak effects. 

In the upper panel of Fig. |l| the normalized top quark 
cross section is displayed in the IS mass scheme 
pTJ versus the cm. energy ^/s, for mj^ — 175 GeV, 
at'^'^\Mz) = 0.118 and Ft = 1.43 GeV. The LL (dot- 
ted lines), NLL (dashed lines), and NNLL (solid lines) 
results are shown for — 0.15, 0.2, and 0.4. From 
a physical point of view the appropriate choice of the 
subtraction parameter v is around as ~ 0.15-0.2. The 
LL curves for f = 0.15,0.2 and 0.4 correspond to the 
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upper, middle and lower lines, respectively. The NLL 
and NNLL curves differ so little on the vertical scale 
of the figure that we refrain from labelling them. At 
NNLL order, the relative variation of the cross sec- 
tion at the peak position is 2% for 0.15 < v < 0.4, 
whereas for -/s = (346,350,351,354) GeV the variation 
is (4%, 2%, 3%, 0.7%). At NLL and LL order the corre- 
sponding variation are 3%, (2%, 3%, 2%, 2%) and 34%, 
(28%, 39%, 27%, 18%), respectively. Overall, the varia- 
tion of the normalization of the cross section for reason- 
able choices of and the shift due to the NNLL or- 
der corrections are an order of magnitude smaller than 
in previous NNLO calculations where threshold masses 
were employed (see Ref. [Q). The improved stability of 
the NNLL cross section is a consequence of the evolu- 
tion of the Wilson coefficients for the potentials and the 
currents. 

In the lower panel of Fig. |^ the contributions to the 
cross section at ^/s ~ 350 Gev coming from (dashed 
line), 5G^ (dotted hue), and the sum of SG^^''''^''' (dot- 
dashed line) are displayed as a function of v. The sum 
of all contributions is represented by the solid line. For 
convenience of presentation we have subtracted the value 



of R 



NNLL 



at V — 0.2 from the dashed and solid curves. 



Whereas the individual contributions vary quite rapidly, 
most notably the Coulomb term G"^ and the l/|k| po- 
tential term SG'', there is a partial cancellation in the 
sum, which varies more slowly. The correlation of the 
Wilson coefficients that leads to this stability is a conse- 
quence of the vNRQCD renormalization group equations 
that account for the correlated soft and ultrasoft run- 
ning and the mixing of Wilson coefficients for v < 1. 
We note that the corrections coming from Vk quickly be- 
come negative for small v, and lead to an instability of 
the NNLL curve at the peak if i' is chosen smaller than 
the Coulombic velocity v ~ 0.15. For v > 0.15 multiple 
energy poles caused by the perturbative treatment of the 
potentials Vs, Vr and Vk and the kinetic energy correction 
do not have to be resummed because the corresponding 
corrections in the binding energies are an order of mag- 
nitude smaller than the top quark width. Ft « 1.5 GeV. 
However, for < 0.15 the residue of the double-pole Vk 
term becomes large and multiple insertions of Vk at the 
peak have to be summed to stabilize the cross section. 
For < 0.1 the value of the Wilson coefficients changes 
rapidly due to the fact that ^jj — mtv^ gets close to 
1 GeV. 

It is instructive to consider the size of some current cor- 
relator corrections from beyond NNLL for v > 0.15. The 
corrections arising from two insertions of Vs are smaller 
than 1% for v > 0.15. In Ref. |l^ the corrections to the 
square of the heavy quarkonium (nS) wave function at 
the origin arising from the emission and rcabsorption of 
an ultrasoft gluon were determined, which can be taken 
as an estimate for ultrasoft corrections to the cross sec- 
tion. To be compatible with our calculations the MS 



parameter used in Ref. IQ has to be replaced by the ul- 
trasoft scale fj,u- For the ground state (n = 1) one finds 
that the corrections amount to about 2% for w . The 
small size of the two corrections just mentioned strength- 
ens our confidence that the 2-3% variation of the nor- 
malization with v at NNLL order represents a realistic 
estimate of the theoretical uncertainties. Using the pre- 
scription given in Ref. |l^] the relative corrections to the 
normalization of the cross section from a 115 GeV stan- 
dard model Higgs boson are 5-8% for energies near the 
threshold. This is larger than the remaining uncertainty 
of the NNLL cross section. 
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